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Quantum transport in disordered magnetic fields is investigated numerically in two-dimensional 
systems. In particular, the case where the mean and the fluctuation of disordered magnetic fields 
are of the same order is considered. It is found that in the limit of weak disorder the conductivity 
exhibits a qualitatively different behavior from that in the conventional random magnetic fields with 
zero mean. The conductivity is estimated by the equation of motion method and by the two-terminal 
Landauer formula. It is demonstrated that the conductance stays on the order of e 2 /h even in the 
weak disorder limit. The present behavior can be interpreted in terms of the Drude formula. The 
Shubnikov-de Haas oscillation is also observed in the weak disorder regime. 



I. INTRODUCTION 



Quantum transport in disordered magnetic fields in 
two dimensions(2D) has been studied widely by numeri- 
cal and by analytical methods. A two-dimensional disor- 
dered system in random magnetic fields with zero mean 
arises in a mean field theory of the fractional quantum 
Hall effect at filling factor v — 1/2 i The magneto trans- 
port around v = 1/2 has then been often analyzed by 
models with random magnetic fields The pos- 
sibility of the Anderson transition in two-dimensional 
system in random magnetic fields has also been studied 
extensively in the last decadei ?'?' 11 ^ 11 ' 12 ^?' 14 Although 
much work has been done to clarify whether or not the 
metallic states exist in such a system which belongs to 
the unitary universality class, the conclusions still re- 
main controversial. Systems in random magnetic fields 
show singularities at the band center in the conductance 
fluctuations 1 ^ and in the density of statesii 4 * 1 ^ These 
singularities are expected to be governed by the chiral 
symmetry of the random magnetic field system. In re- 
cent experiments, two-dimensional electron systems in 
random magnetic fields have been constructed and the 
magneto transport in such systems has been measured^ 
The random magnetic fields with zero mean are realized 
by attaching small magnets on the layer parallel to the 
2D electron gas in a semiconductor heterostructure. It 
is then found that the magneto resistance exhibits simi- 
lar structure to that observed in the fractional quantum 
Hall effect around v = 1/2. It is thus important to un- 
derstand the transport in random magnetic fields both 
theoretically and experimentally. 

In the present paper, we focus on another aspect of the 
transport in disordered magnetic fields. In the previous 
paperji& we have investigated the magneto transport in 
random magnetic fields and found that the conductivity 
is insensitive to the magnitude of the fluctuation of the 
random magnetic fields if the mean value B of random 
magnetic fields is set to be of the same order of its fluc- 



tuation SB. This is in contrast to the case of B = 
where the conductivity diverges as oc (SB)" 2 in the limit 
of SB — s- 0. It is thus an interesting problem to clarify 
whether the system exhibits this insensitivity even in the 
limit of zero magnetic fields. We thus focus on the trans- 
port properties in the case where SB w B, particularly 
in the limit of weak magnetic fields. The characteris- 
tic feature of such systems is that the magnetic fields 
have almost the same sign. We therefore consider, as a 
simplified model, the two-dimensional system in random 
magnetic fields with a fixed sign in order to analyze the 
case SB k, B. 

Generally, it is expected that the scattering mecha- 
nism in random magnetic fields with a fixed sign is qual- 
itatively different from that in the conventional random 
magnetic field with zero mean. The system has, for ex- 
ample, no snake state near the zero magnetic field lines, 
which plays an important role in the semi-classical the- 
ory of the transport in the conventional random magnetic 
fieldsAi It would thus be useful to study the present sys- 
tem also for the further understanding of the quantum 
transport in the conventional random magnetic fields. 
On the other hand, in the strong magnetic field limit 
B ^> SB it has been argued and demonstrated that the 
statistical properties belong to the same universality class 
as the quantum Hall transition in two dimensions^iSi 
The float of extended states toward the limit of the weak 
magnetic field B <C SB has then become an impor- 
tant issued 0-21 Clarifying the transport properties of the 
present case (B ss SB) , particularly in the limit of weak 
magnetic fields, would be important also in considering 
the connection between these two limiting cases. 

We evaluate numerically the conductance in random 
magnetic fields with a fixed sign by the equation of mo- 
tion method and by the Landauer formula. In the equa- 
tion of motion method, the conductance is obtained by 
examining the electron diffusion directly. This method 
has an advantage that it applies to very large systems. 
On the other hand, in the case of the Landauer formula, 
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large numbers of samples can be considered, although 
system-sizes are limited. With these two independent 
numerical method, we calculate the longitudinal conduc- 
tivity and find that the conductance exhibits a qualita- 
tively different behavior from that in the conventional 
random magnetic fields, particularly, in the limit of weak 
disorder. We also show that the present observations are 
not changed significantly by introducing a small number 
of magnetic fields with the opposite sign. 

II. MODEL 

We adopt a model described by the following Hamil- 
tonian on the square lattice 

H= J2 VeMMi,i)CjCj, (1) 

<i,j> 

where C\(Ci) denotes the creation(annihilation) opera- 
tor of an electron on the site i. The phases {Oij} are 
related to the disordered magnetic fluxes {</>j} through 
the plaquette (i, i + x, i + x + y, i + y) as 

&i,i+x +9i+x,i+x+y +9i+x+y,i+y +&i+y,i = ~2TT(f>i J (pQ (2) 

where 0o — h/\e\ stands for the unit flux and x(y) de- 
notes the unit vector in the a;(y)-direction. The fluxes 
{4>i \ are assumed to be distributed independently in each 
plaquette. The probability distribution P(<fi) of the flux 
cj> is given by 

pu) = { y^ f { °i ° ^ ^ ftrf . (3) 

y ' 10 otherwise 

The mean and the variance of the distribution are ac- 
cordingly given by 

(4>i) = -^<t>Q and (Mi)-(4>i)(4>j) = ^<ft8 i j, (4) 
respectively. 

III. ELECTRON DIFFUSION 

In order to verify that the above system has a diffu- 
sive regime, we first examine the diffusion of electrons by 
means of the equation of motion method. For this, we nu- 
merically solve the timc-dcpcndent Schrodingcr equation 
by making use of the decomposition formula for exponen- 
tial operators^ The basic formula used in the present 
paper is the forth order formula 

exp(a:L4i+- • ■+A n ]) = S(xpf S(x{l-4p))S(xp) 2 +0(x 5 ), 

(5) 

where 

S( X ) = e xM ' 2 ■ ■ ■ e xA n - 1 /2 e xA ne xA n _ 1 /2 . . . e x Al /2^ (g) 



The parameter p is given by p — (4 — 4 1 / 3 ) -1 and 
Ai,...,A n are arbitrary operators. We divide the Hamil- 
tonian into five parts as in the previous papers^ 3 - so that 
each part is represented as the product of 2 x 2 subma- 
trices. By applying this formula to the time evolution 
operator U(t) = exp(—iHt/Ti), we obtain 

U(St) = U 2 (pSt) 2 U 2 ((l - 4 P )St)U 2 (pSt) 2 + 0(6t 5 ) (7) 

with 

U 2 ( X ) = e xH ' /2 ■ .. e ^/2 e xH, e x Hl /2 . .. e x^/ 2j (g) 

where H = Hi + • • • + H$. It is to be noted that the U 2 
can be expressed in an analytical form while the original 
evolution operator U can not be evaluated exactly with- 
out performing the exact diagonalization of the whole 
system. This method has already been applied to the 
various cases of 2D 18,23,24 and 3D 25 disordered systems. 

Numerical calculations are performed in the system 
of the size 999 x 999 with the fixed boundary condi- 
tion. All the length scales are measured in units of the 
lattice constant a. To prepare the initial wave packet 
with energy E, we numerically diagonalize the subsys- 
tem (21 x 21) located at the center of the whole system 
and take as the initial wave function the eigenstate whose 
energy is closest to E. The single time step St is set 
to be 5t = 0.2h/V and at least five realizations of ran- 
dom magnetic fields are considered. With this time step, 
the fluctuations of the expectation value of the Hamilto- 
nian is safely neglected throughout the present simulation 
(t < 2000h/V). We observe the second moment defined 
by 

(r 2 (t)) c ^(r 2 (t))-(r(t)) 2 (9) 

with 

{r n {t)) = J r'Mft^-MrlV'M)! 2 (10) 

where ^(r, t) denotes the wave function at time t. In the 
diffusive regime, the second moment is expected to grow 
in proportion to t 

{r 2 ) c = 2dDt, (11) 

where the diffusion coefficient is denoted by D and d is 
the dimensionality of the system. The diffusion coeffi- 
cient D is related to the conductivity by the Einstein 
relation a = e 2 Dp. Here p denotes the density of states. 
It is estimated by the Green function method 26 for strips 
with the width up to 30. 

The second moment (r 2 ) c for h T { = 0.04 as a function 
of time t is shown in Fig. 1. The fermi energy is as- 
sumed to be E/V = 0.5. It is clearly seen that, above 
a certain length scale, the second moment grows in pro- 
portion to the time t. In this regime, it is natural to 
expect that the system is diffusive. The length scale, 
above which the diffusive behavior is observed, must be 
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FIG. 1: The second moment averaged over 10 realization of 
random magnetic fields as a function of time t for E/V = 0.5 
and h v { — 0.04. The fluctuation around the mean value is 
plotted at t/(n/V) = 500, 1000, 1500 and 2000. 



FIG. 2: The conductivity as a function of /i r f for E/V = 0.5. 
For h r { < 0.1, the conductivity stays on the order of the 
conductance quantum e 2 /h. 



related to the typical cyclotron radius of the present sys- 
tem. Since the fermi energy is fixed close to the band 
center, it would vary almost inversely proportional to h r { . 
In fact, the length scales are estimated as ~ 38, ~ 49, and 
~ 70 for hrt = 0.04, 0.03 and 0.02, respectively. Below 
these length scales, the system is in a ballistic regime. In 
estimating the diffusion coefficient, we discard the data 
in the ballistic regime. It should be noted that, in the 
present regime of h T {, the effective mean free path due 
to the fluctuation of magnetic fields is much larger than 
the cyclotron radius determined by the mean value of the 
magnetic fields. For instance, for the case of h T { = 0.02, 
the effective mean free path is estimated to be around 
2100api& whereas the cyclotron radius is 40a. 

The estimated conductivity is plotted as a function of 
h T { in Fig. 2. It is rather remarkable that the conductiv- 
ity stays on the order of the conductance quantum e 2 /h 
even for the small values of h T {. This feature is in contrast 
to the case of the conventional random magnetic fields, 
where the conductivity is likely to diverge inversely pro- 
portional to the square of the fluctuation of the random 
magnetic fields. These facts mean that, provided that 
the sign of the magnetic fields is fixed, the conductivity 
is insensitive to the strength of its fluctuation, and im- 
plies that the zero field limit h T { —> does not coincide 
with the hrf = case and thus is a singular limit. 

IV. CONDUCTANCE 

In order to investigate this small conductivity in the 
weak fluctuation regime, we examine also the conduc- 
tance in the two-terminal geometry based on the Lan- 
dauer formula 

G = —Tr(TT^), (12) 
h 

where T denotes the transmission matrix through the 
system. We consider the L by L square system with 



two leads connected to both sides of the system. No 
magnetic field is assumed in the leads. For this system, 
we adopt the transfer matrix formalism 2 ^ for evaluating 
the transmission coefficient numerically. The maximum 
system-size considered in the present work is 160 by 160 
and the average over 1000 realizations of random field 
configurations is performed. Since the mean value of the 
magnetic fields is not zero in the present system, the ef- 
fect of the edge states is, in general, important for the 
transport properties. We therefore estimate the conduc- 
tance both with the fixed boundary conditions and with 
periodic boundary conditions in the transverse direction. 

In Fig. 3, the conductances for the fixed boundary con- 
ditions for various system-sizes are shown as a function 
of h r {. It is clearly seen that the conductance is likely to 
diverge as h T { — > 0, which is qualitatively different from 
the behavior obtained by the analysis of the electron dif- 
fusion. It is also seen that the value of the conductivity 
is much larger than that estimated from the electron dif- 
fusion. It is expected that these differences come from 
the fact that the edge states, which are absent in the 
simulation of the electron diffusion, exist in the case of 
the fixed boundary conditions. The observation that the 
conductance is almost inversely proportional to hrf (Fig. 
3) suggests that it is the manifestation of the fact that 
the number of edge states increases as fiZ . 

In order to remove the effect of edge states, we evaluate 
the conductivity under the periodic boundary conditions 
in the transverse direction (Fig. 4). It is seen that, in 
contrast to the case of the fixed boundary conditions, 
the conductance stays on the order of the conductance 
quantum even for the small values of h T f. It should be 
emphasized that for small h T { the conductance decreases 
as the size of the system is increased. We thus recover 
the results obtained from the electron diffusion. This 
suggests that it is essential to remove the effect of edge 
states for observing the small conductance in the random 
magnetic fields with a fixed sign. 

In the case of the periodic boundary conditions, the 
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FIG. 3: The conductance under the fixed boundary conditions 
for E/V = 0.5 and L — 60(open circles), 80(open squares), 
100(open triangles), 120(solid triangles), 140(solid squares) 
and 160(solid circles). The average over 1000 realizations of 
random fields is performed. 
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FIG. 4: The conductance under the periodic boundary con- 
ditions for E/V = 0.5 and L = 60(open circles), 80(open 
squares), 100(open triangles), 120(solid triangles), 140(solid 
squares) and 160(solid circles). The average over 1000 real- 
izations of random fields is performed. 



total flux through the system is adjusted to be an integer 
multiple of the flux quantum h/\e\ in order to ensure that 
the flux through the leads to be zero. Practically, we 
divide the excess flux into small pieces and subtract these 
pieces randomly from the fluxes already distributed in 
the system as long as they are positive, so that the above 
condition is satisfied. This procedure produces a weak 
correlation in random fluxes which should be negligible 
for large system sizes. 

More detailed calculations of the conductance for the 
periodic boundary conditions are shown in Fig. 5. It is 
clearly seen that the conductance oscillates periodically 
as a function of l//i r f and its period is independent of 
the system-sizes. From our numerical data (Fig. 5), the 
period in l//i r f is estimated to be A(l//i r f) = 1.33. The 
same type of oscillation is observed also for E/V = 0.3 
and for E/V = 0.7 with different periods. The periods 
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FIG. 5: The conductance under the periodic boundary con- 
ditions for E/V = 0.5 and L — 80(open squares), 100(open 
triangles), 120(solid triangles) as a function of l/h T f. The 
period in l//i r f is estimated to be 1.33. 

in l/h r [ are estimated to be 1.18 and 1.41, respectively. 
Here we show that these oscillations can be understood 
as the Shubnikov-de Haas effect^ In general, the period 
of the Shubnikov-de Haas effect for tight-binding lattice 
models is given by 

4tt 2 

A(l/0ext) = ~2q, 0oxt = ^cxt/00 (13) 

where </> GXt denotes the external flux through the pla- 
quette and S the area of the fermi surface of the two- 
dimensional electron system. By evaluating S in the sys- 
tem without magnetic fields, we obtain A(l/0 cxt ) = 2.36, 
2.59 and 2.83 for E/V = 0.3, 0.5 and 0.7, respectively. 
Since, in the present system, the mean value of the 
magnetic field (4>/4>o) is given by h T {/2, we can iden- 
tify (f) ext to be h r f/2. We then find the relation that 
A(l/</> e xt) = 2A(l//i r f). With this relation, it is clear 
that the periods of the Shubnikov-de Haas effect eval- 
uated by ea. l|13|) are in excellent agreement with those 
of the oscillation observed in our numerical data. It is 
rather remarkable that the present system exhibits such 
a clear Shubnikov-de Haas effect even though the flux 
distributes uniformly from to h T f. A smaller oscillation 
than that of the conductance is observed in the density 
of states (Fig. 6), which seems to be consistent with the 
analytical results&22i for the conventional Shubnikov-de 
Haas oscillation of 2D systems in magnetic fields. 

Apart from this oscillation, the conductance shows a 
smooth cross over from the small h r { regime to the htf — 1 
limit. It is to be noted that due to the periodicity of 
the phases 9% t j in the hamitonian, the h r { = 1 limit is 
identical to the case of the conventional random magnetic 
field with zero mean. Our numerical data suggest that 
this cross over takes place around h r i = 0.1 ~ 0.2. With 
this cross over, the conductance is increased to a value 
~ 1.5e 2 /h (Fig. 5) as h T { is increased. This enhancement 
of conductance in the regime h T { > 0.1 is also observed in 
the the conductivity estimated by the equation of motion 
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FIG. 6: The density of states for E/V = 0.5. 



method (Fig. 2). 



FIG. 7: The conductance for E/V = 0.5 in the case 
where the random fluxes <fi/<j>o are distributed in the range 
[— 0.05/i r f , 0.95/irf]- Open squares, open triangles, solid trian- 
gles represent L = 80, 100 and 120, respectively. 



V. DISCUSSION 

We have investigated the longitudinal conductivity in 
the two-dimensional disordered magnetic fields with a 
fixed sign. It has been demonstrated that the system 
shows a slow diffusion and has accordingly a small con- 
ductance which is on the order of e 2 /h and is insensitive 
to the magnitude h T { of the random fields. It would be 
useful to interpret these results in terms of the Drudc 
formula a = 00/(1 + uj c t2 )i where w c and r denote the 
cyclotron frequency and the relaxation time, respectively. 
The conductivity for oj c = is denoted by tjQ. Here it 
should be kept in mind that the present system is equiv- 
alent to the system in the conventional random magnetic 
field distributed from —h r {/2 to htf/2, plus an additional 
uniform magnetic field h T {/2. Namely, we can regard it 
as the conventional random field system in an external 
uniform magnetic field. From this point of view, we as- 
sume that the relaxation time r is determined simply 
by the scattering due to the conventional random fields 
with zero mean and is insensitive to the additional uni- 
form magnetic field. We also assume that the effect of 
the additional uniform magnetic field appears only in tu c . 
It is then deduced^ that ctq oc t oc h~ 2 and lo c oc h r f, 
and therefore that for h T { <C 1, 



a/(e 2 /h) 



<?o/(e 2 /h) 
(1+^ T 2) ~ (1 + 



A iKf 2 



(14) 



where A% and A2 are constants independent of /i r f. If we 
take the limit as h Y f — > 0, we obtain that a/(e 2 /h) — > 
A\jAi. The parameters A\ and A2 are basically deter- 
mined by the fermi energy E&L For instance, A\ has 
been estimated to be 0.95 for E/V = 0.5 M The present 
argument based on the Drude formula, though it does 
not account for the Shubnikov-de Haas effect, seems to 
account for the fact that for h T { -C 1 the conductance 
stays on the order of e 2 /h and is insensitive to h T {. In 



this context, the feature specific to the present system is 
that both t and (toj c ) 2 scales in the same way as h~ 2 . 

Finally, we emphasize that these present transport 
properties are not restricted to the case of the strictly 
positive random magnetic fields. These properties can 
be observed in the case where the random magnetic 
fields are almost positive. In Fig. 7, the conductance 
is shown for the case where the random fluxes {(f>i/<f>o} 
are distributed in the range [— 0.05/i r f, 0.95/i r f]. It is clear 
that the Shubnikov-de Haas effect is also observed even 
though the magnetic fields are not strictly positive. This 
suggests that these properties are general features of sys- 
tems in the random magnetic fields whose mean and fluc- 
tuation are of the same order. 

In summary, we have studied numerically the transport 
properties in two-dimensional random magnetic fields 
with a fixed sign. In particular, the conductivity in the 
limit of small fluctuation has been investigated. In the 
simulation of the diffusion of electron wave functions, 
it is observed that the system shows diffusive behavior 
in larger length-scales than the typical cyclotron radius. 
The diffusion turns out to be very slow and the conduc- 
tivity evaluated from the diffusion coefficient stays on 
the order of the conductance quantum insensitive to the 
magnitude of fluctuation of the magnetic fields. This 
behavior is also observed in the two-terminal conduc- 
tance evaluated by means of the Landauer formula under 
the periodic boundary conditions. A clear Shubnikov-de 
Haas effect has been observed for a wide range of h T {. 
For a system with edge states, the contribution to the 
conductance from the edge channels is much larger than 
that from the bulk, and hence these properties are not 
observable. In this sense, these peculiar features are the 
bulk properties of the positive random magnetic fields. 
It has been argued that the singularity in the limit of 
/i r f — ► can be understood within the framework of the 
Drude formula. 
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